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ABSTRACT 

Mechanical  behavior  of  granular  media  is  important  in  many  fields  of 
studies  such  as  soil  mechanics,  powder  mechanics,  and  ceramic  mechanics. 

The  mechanical  behavior  of  granular  media  has  been  studied  by  borrowing  the 
stress-strain  models,  such  as  elastic,  elasto-plastic,  or  plastic  models, 
developed  for  continuum  materials .  These  continuum  models  consider  neither 
the  discrete  nature  nor  the  deformation  mechanism  of  granular  materials. 

A  more  rational  approach  should  be  one  that  considers  the  granular 
system  as  an  assemblage  of  particles.  The  stress-strain  behavior  for  a 
granular  material  is  define  for  a  representative-volume  which  consists  of  a 
sufficiently  large  number  of  particles  to  be  representative  of  the  material. 
When  subjected  to  loading,  the  deformation  of  the  granular  system  results 
from  particle  deformation  as  well  as  slip  between  particles. 

Along  this  line,  the  aim  of  this  research  is  to  investigate  the 
constitutive  behavior  of  granular  assemblies  using  a  micro-mechanics 
approach,  taking  into  account  the  non-uniform  strain  within  the 
representative-volume.  The  nature  of  this  investigation  is  focused  on 
theoretical  development.  Formulation  of  macroscopic  constitutive  theory, 
firmly  founded  in  micromechanics  of  particle  interaction,  provides  an 
improved  understanding  and  description  of  granular  deformation  behavior. 

The  developed  theory  is  evaluated  by  results  obtained  from  computer 
simulation  and  experiments. 
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CHAPTER  1 

SOJMARY  OF  THE  PROJECT 


1.1  RESEARCH  OBJECTIVES 

The  overall  otjective  of  this  research  is  to  develop  a  constitutive 
model  for  granular  materials  under  finite  strains  with  consideration  of 
particle  sliding  and  the  consequent  fabric  change.  The  specific  objective  of 
this  research  is  focused  on  the  development  of  a  stress-strain  theory  that 
accounts  for  the  microstructure  and  the  non-uniform  strain  field  of  the 
particle  assembly.  The  developed  theoretical  model  is  evaluated  by  results 
obtained  from  computer  simulation  and  experimental  tests. 

1.2  ACCOMPLISHMENTS 

A  general  stress-strain  theory  is  developed  for  granular  soil  which 
explicitly  accounts  for  the  microstructure  and  the  micro-deformation 
mechanisms  of  soil,  such  as  sliding  and  separation  of  particles.  The  theory 
adopts  a  statistical  mechanics  approach  in  which  the  geometric  structure  of 
the  particulate  material  is  characteristically  represented  by  the 
configurations  of  a  set  of  micro-elements.  Each  micro-element  is  defined  as 
a  particle  group  consisting  of  a  center  particle  and  its  neighbor  particles. 
Following  the  concepts  of  statistical  mechanics,  the  stress-strain  behavior 
of  the  material  can  be  interpreted  as  the  average  behavior,  in  a 
probabilistic  sense,  of  the  set  of  micro-elements. 

Most  of  the  research  results  have  been  published  in  Journals  and 
Conference  proceedings.  A  brief  summary  of  the  results  is  given  here.  The 
reference  number  in  parenthesis  noted  in  this  summary  is  referring  to  the 


sequential  number  in  the  publication  list  given  in  section  1.3.3  (i.e., 

Publications  under  this  grant) . 

1.2.1  Theoretical  development 

a)  Continuum  concepts  of  stress  and  strain  were  developed  for  a  micro¬ 
element  in  terms  of  the  discrete  variables,  such  as  contact  forces  and 
particle  movements  (2,  4,9). 

b)  A  new  definition  of  strain  was  derived  from  kinematics  of  particles. 
Unlike  that  in  classic  mechanics,  the  strain  for  granular  material  is 
asymmetric.  The  symmetric  part  of  the  strain  is  the  usual  Cauchy  strain. 
The  skew-symmetric  part  of  the  strain  represents  the  spin  of  particles 
(7,  4,  9). 

c)  Constitutive  theory  was  developed  considering  the  asymmetric  stress 
tensor  resulted  from  rolling  resistance  between  particles.  The  granular 
material  is  derived  to  be  a  more  generalized  continuum  in  which  couple 
stress  can  transmit.  The  effects  of  rolling  resistance  between  particles 
on  the  deformation  behavior  of  material  are  discussed  (9,  12,  14) . 

d)  A  general  stress-strain  theory  for  random  packing  of  multi-sized  granules 
has  been  developed.  The  stress-strain  behavior  is  based  on  the 
properties  at  contacts  between  particles  (3,  9,  16) . 

e)  The  theory  has  been  extended  to  large  strain  conditions  including  the 
mechanisms  of  sliding  and  separation  of  particles  associated  with  non- 
uniform  strain  field  (17,  18) . 


1.2.2  Computer  Simulation 
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a)  The  theory  has  been  verified  by  computer  simulation  for  assemblies  of 
spheres  and  disks  (5,  11,  14),  and  by  experimental  results  from  rod 
assemblies  tested  in  a  direct  shear  box  (2,  3,  7) . 

b)  The  theory  has  been  applied  to  boundary  value  problems  using  an 
innovative  numerical  procedure  similar  to  finite  element  formulation  (12, 
13,  14)  . 

1.2.3  Experimental  Evaluation 

a)  The  theory  has  been  applied  to  study  the  behavior  of  sand  (6,  15,  18)  and 
cemented  sand  (10)  on  initial  moduli,  secant  moduli  and  damping  ratio, 
compared  with  experimental  results  for  low  amplitude  loading  conditions. 

b)  The  theory  has  been  studied  for  the  effects  of  fabric  (packing  structure) 
on  the  class  of  material  symmetry  (2,  8,  18) . 

c)  The  extended  theory  for  large  strain  conditions  including  the  mechanisms 
of  sliding  and  separation  of  particles  (17,  18)  is  evaluated  by  comparing 
the  predicted  stress-strain  behavior  with  that  measured  from  experiments 
on  different  types  of  sands  under  various  loading  conditions  (18) . 

In  summary,  the  comparisons  between  the  predicted  and  measured  results 

for  small  and  large  strain  conditions  show  that  the  model  has  the  foilcwing 

capabilities : 

a)  The  model  is  capable  of  predicting  macroscopic  behavior  such  as  moduli 
and  failure  of  granular  material  based  on  the  microscopic  contact 
properties . 

b)  The  effect  of  microstructure,  such  as  inherent  anisotropy,  can  be 


explicitly  accounted. 


c)  The  model  accounts  for  the  behavior  associated  with  particle  interaction 
and  sliding  at  contact,  thus  capable  of  predicting  behavior  such  as 
stress-induced  anisotropy,  path  dependency,  locked-in  stress,  plastic 
flow,  dilatancy,  friction  loss  (damping) ,  and  non-coaxial  behavior  under 
rotation  of  principal  stress. 

d)  The  model  is  applicable  to  materials  with  cemented  inter-particle 
properties . 
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CHAPTER  2 

THEORETICAL  DEVELOPMENT 


The  intrinsic  feature  of  granular  materials  is  their  particulate 
structure.  The  micro-scale  interaction  between  particles  at  contacts 
controls  the  macroscopic  response.  Mechanical  behavior  of  granular 
material  can  be  analyzed  by  two  approaches  which  account  for  the 
particulate  structure  and  the  particle  interactions,  namely,  the 
computer  simulation  approach  and  the  micromechanics  approach.  The 
computer  simulation  approach  solves  for  the  deformation  of  an  assembly 
based  on  the  governing  equations  for  the  movement  of  each  particle 
interacting  with  its  surrounding  particles.  Along  this  line  of 
approach,  methodologies  can  be  found  in  the  work  by  Serrano  and 
Rodriguez-Ortiz  (1973)  and  Kishino  (1988)  for  quasi-static  condition  and 
in  the  work  by  Cundall  and  Strack  (1979)  for  dynamic  condition. 

However,  this  approach  is  cumbersome  for  systems  composed  of  large 
number  of  particles.  For  example,  the  number  of  particles  is  in  the 
order  of  hundreds  of  thousands  for  a  cubic  inch  of  sand.  Hence  it  is 
desirable  to  represent  the  discrete  system  with  a  more  tractable 
continuum  model. 

In  the  micromechanics  approach,  the  deformation  behavior  of  the 
assembly  is  described  by  the  continuum  concepts  of  stress  and  strain. 

The  constitutive  relationship  for  a  granular  assembly  is  obtained  by 
considering  the  material  structure  and  particle  interaction.  A  number 
of  studies  have  been  attempted  along  this  line  of  approach.  For 
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example,  work  can  be  found  in  Duffy  (1959),  Duffy  and  Mindlin  (1957), 
Deresiewicz  (1958),  and  Makhlouf  and  Stewart  (1967)  for  regular 
packings,  and  Digby  (1981),  Walton  (1987),  Jenkins  (1987),  Chang  (1988), 
and  Bathurst  and  Rothenburg  (1988)  for  random  packings.  However,  all 
the  aforementioned  models  are  limited  to  the  assumption  that  there  is  no 
sliding  at  contacts  during  deformation,  thus  only  applicable  to  the 
conditions  of  small  deformation. 

In  this  study,  we  endeavor  to  extend  the  constitutive  relationship 
for  granular  systems  taking  into  account  the  effect  of  particle 
separation  and  sliding  under  large  deformation.  One  of  the  major 
obstacles  to  account  for  particle  sliding  under  large  deformation  is  the 
non-uniform  strain  field  (Chang  1989,  Chang  and  Misra  1990b) .  In  this 
chapter,  we  tackle  this  problem  by  introducing  a  distributive  law  which 
describes  the  heterogeneous  field  of  strain. 

In  what  follows,  we  first  briefly  review  the  previous  work  on 
granular  mechanics.  Then  we  describe  the  current  approach  of  treating 
the  material  at  three  levels,  namely,  representative  unit,  micro-element 
and  inter-particle  contact.  On  this  basis,  the  constitutive  laws  for 
each  level  are  derived  and  the  overall  stress-strain  relationship  is 
described  in  terms  of  inter-particle  contact  behavior.  The  performance 
of  this  model  compared  with  experimental  results  are  discussed  in 
chapters  3  and  4 . 


2.1  REVIEW  OF  PREVIOUS  WORK  ON  GRANULAR  MECHANICS 
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We  categorize  the  work  on  constitutive  mechanics  of  granular 
material  into  two  areas,  namely,  1)  mobilized  plane  concept,  and  2) 
micromechanical  concept. 

2.1.1  Mobilized  plane  concept 

The  concept  of  mobilized  plane  was  first  introduced  by  Coulomb 
(1773)  to  describe  the  resistance  due  to  internal  friction  between 
particles.  Based  on  this  concept,  Rankine  (1057)  developed  theory  of 
active  and  passive  pressures  and  Caquot  (1934)  derived  a  solution  of 
limiting  equilibrium  for  an  ideal  granular  wedge. 

Taylor  (1948)  and  Bishop  (1950)  differentiated  the  sliding  within 
granular  material  from  that  of  two  solid  blocks  by  addressing  the 
mechanism  of  interlocking  between  particles  of  granular  material. 
Interlocking  of  particles  restricts  the  degree  of  mobilization  thus 
attributes  to  a  large  portion  of  shear  strength.  As  a  result  of  the 
interlocking,  dilatancy  occurs  with  sliding. 

Referring  to  the  notion  of  dilatancy,  Newland  and  Allely  (1957) 
suggested  that  the  relative  direction  of  sliding  between  two  blocks  is 
in  general  not  parallel  to  the  mobilized  plane,  rather  it  is  inclined  at 
an  angle.  This  mechanism  leads  to  the  phenomenon  of  dilatancy  and, 
subsequently,  to  the  difference  between  peak  and  residual  stresses. 

Rowe  (1962)  postulated  an  energy  ratio  criterion  to  derive  the  relative 
direction  of  sliding  between  two  blocks  of  particles  in  a  random 
assembly  subjected  to  a  triaxial  loading  condition.  Based  on  this 
derivation,  the  angle  between  sliding  direction  and  mobilized  plane  is  a 
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function  of  applied  stress.  Horne  (1965,  1969)  studied  the  Rowe's 
energy  postulate  in  a  more  rigorous  way  by  considering  the  sliding 
between  pairs  of  particles  in  a  random  assembly.  It  is  shown  that  the 
average  of  these  sliding  directions  can  be  obtained  by  maximizing  the 
energy  transmission  ratio. 

Along  this  line,  by  accounting  not  only  the  directions  but  also  the 
magnitudes  of  sliding  over  mobilized  planes,  general  stress-strain 
relationship  have  been  derived  by  many  investigators,  for  example, 
Matsuoka  (1974),  de  Josselin  de  Jong  (1977),  Prat  and  Bazant  (1991)  and 
Chang  et  al.  (1989c) . 

In  the  mobilized  plane  concept,  it  is  necessary  to  describe 
phenomenologically  the  averaged  behavior  of  particles  on  a  mobilized 
plane.  The  microscopic  structure  and  particle  interaction  are  not 
directly  considered.  It  is  therefore  desirable  to  chose  an  alternative 
approach  based  on  micromechanicai  concept  which,  at  a  more  fundamental 
level,  considers  explicitly  the  laws  of  mechanics  governing  the 
interaction  of  particles  at  contacts.  The  two  concepts  are 
schematically  illustrated  in  Fig.  2.1. 

2.1.2  Micromechanicai  concept 

The  earliest  attempt  based  on  this  micromechanics  approach  was  by 
Duffy  and  Mindlin  (1957) .  They  presented  a  stress-strain  relation 
derived  for  a  medium  composed  of  a  face-centered  cubic  array  of  elastic 
spheres  in  contact.  Following  their  approach,  stress-strain 
relationship  have  been  developed  by  Deresiewicz  (1958)  for  simple  cubic 


Sliding  Plane 


A.  Mobilized  Plane  Concept 


B.  Micromechanical  Concept 
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Fig.  2.1  Two  concepts  in  constitutive  mechanics  of  granular  material:  a) 
mobilized  plane  concept,  and  b)  micromechanical  concept. 
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array,  and  by  Makhicuf  and  Stewart  (1967)  for  cubical-tetrahedral  and 
tetragonal  sphenoidal  arrays. 

Failure  of  a  face-centered  packing  has  been  investigated  by 
Thurston  and  Deresiewicz  (1959)  and  Rennie  (1959) .  Rennie  applied  this 
results  to  the  prediction  of  failure  of  random  arrays,  arguing  that  the 
random  assembly  consists  of  many  face-centered  miniature  arrays  in 
various  orientations.  Failure  initiates  when  the  most  unfavorably 
oriented  miniature  arrays  fail,  causing  the  failure  of  the  entire 
assembly  of  the  particles.  This  concept  of  failure  is  however 
contradicted  by  experimental  evidence. 

More  recently,  with  the  aid  of  the  knowledge  gained  in  mechanics 
for  composite  material,  the  micromechanics  approach  for  granular 
material  has  rapidly  developed  and  shown  fruitful  results. 

Christoffersen  et  al.  (1981)  represented  macroscopic  stress  of  an 
assembly,  on  the  basis  of  the  principle  of  virtual  work,  in  terms  of 
microscopic  contact  forces.  Stress-strain  models  for  packings  of  equal 
sized  spheres  under  axisymmetric  loading  conditions  have  been  developed 
by  Walton  (1987)  for  isotropic  packings  and  by  Jenkins  (1987)  for 
idealized  anisotropic  packings.  A  general  stress-strain  theory  for 
random  packing  of  multi-sized  granules  has  been  developed  by  Chang 
(1988)  and  Chang  et  al.  (1989a) .  The  theory  has  been  verified  by 
computer  simulation  of  disks  (Chang  and  Misra  1989b)  and  by  experimental 
results  from  rod  assemblies  (Chang  and  Misra  1989a) .  The  theory  has 
been  applied  to  study  the  behavior  of  sand  (Chang  et  al.  1969b)  and 
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cemented  sand  (Chang  et  al.  199Q) ,  with  discussions  of  the  fabric 
effects  (packing  structure)  on  initial  moduli  (Chang  and  Misra  1990a) . 

It  is  noted  that  all  the  aforementioned  stress-strain  models  for 
randomly  packed  granular  material  are  limited  to  small  strain  conditions 
without  the  consideration  of  particle  sliding.  In  the  following,  we 
aim  to  extend  the  previous  model  taking  into  the  consideration  of 
mechanisms  of  sliding  and  separation  between  particles. 

2.2  THREE-LEVEL  MICROMECHANICAL  APPROACH 

Granular  material  can  be  viewed  at  three  levels,  namely,  1) 
representative  unit,  2)  micro-element,  and  3)  contact.  A  representative 
unit  is  defined  as  an  assembly  which  contains  large  number  of  particles 
to  be  representative  of  the  granular  materials.  Each  particle  in  the 
representative  unit  is  in  contact  with  several  neighboring  particles. 
These  particles  form  a  particle  group,  termed  as  micro-element.  A 
micro-element,  corresponding  to  each  particle,  is  an  elementary  unit  at 
microscopic  level.  There  are  several  contacts  between  the  center 
particle  and  its  neighbors.  A  contact  between  a  pair  of  particles  is 
regarded  as  the  basic  unit  of  granular  material.  A  schematic 
representation  of  the  three  levels  of  granular  material  is  shown  in  Fig. 
2.2. 

Figure  2.3  shows  the  concept  of  micromechanics  approach  used  in  this 
chapter.  In  such  approach,  the  constitutive  behavior  are  defined  for 


three  levels: 


MICRO-ELEMENT 


REPRESENTATIVE  UNIT 


Fig.  2.2  Schematic  representation  of  three  levels  of  granular  material. 
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Fig.  2.3  Micromechanics  approach  for  modelling  the  mechanical  behavior  of 


granular  material. 
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1)  At  contact  level,  the  constitutive  law  relates  contact  force  and 
relative  movement  between  two  particles.  At  this  level,  continuum 
concept  has  not  yet  been  introduced.  We  are  dealing  with  the 
translation  and  rotation  of  discrete  particles  which  result  to  the 
micro-scale  discontinuity  due  to  mechanisms  such  as  sliding  and 
separation. 

2)  At  micro-element  level,  the  constitutive  law  relates  the  stress 
and  strain  for  the  micro-element.  The  continuum  concepts  of  stress  and 
strain  are  now  introduced.  We  seek  to  derive  the  stress-strain 
relationship  for  a  micro-element  based  on  contact  behavior.  To 
accomplish  this  objective,  it  is  necessary  to  establish:  a)  the 
relationship  between  stress  and  contact  forces,  and  b)  the  relationship 
between  particle  movement  and  strain. 

3)  At  representative  unit  level,  constitutive  law  relates  the 
overall  stress  and  overall  strain.  The  behavior  of  the  representative 
unit  is  obtained  by  averaging  the  behavior  of  micro-elements  utilizing 
the  principle  of  volume  average.  A  concept  of  'distributive  law  of 
strain'  is  introduced  to  facilitate  the  averaging  process.  Therefore, 
the  macro-behavior  of  an  assembly  can  ultimately  be  derived  from  the 
micro-behavior  of  a  contact. 

In  what  follows,  the  theoretical  development  of  constitutive 
behavior  are  described  for  the  three  levels:  1)  contact,  2)  micro¬ 
element,  and  3)  representative  unit.  The  objective  is  to  derive  the 
constitutive  law  of  a  representative  unit  based  on  the  behavior  of 
inter-particle  contact. 
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2.3  BEHAVIOR  AT  CONTACT  LEVEL 

For  materials  such  as  granular  soil,  powder,  and  ceramic,  granules 
are  relatively  rigid  and  the  deformation  of  granules  occurs  mostly  at 
contacts.  Thus  the  granular  material  is  envisioned  to  be  rigid 
particles  connected  by  imaginary  elastic-plastic  springs  which  allow 
sliding  and  separation  of  particles. 

There  are  two  modes  of  movement  for  a  particle:  translation,  u., 
and  rotation,  on  .  Based  on  the  kinematics  of  two  rigid  particles  of 
convex  shape,  the  relative  displacement  A6™  and  relative  rotation  A0™ 
between  particle  'n'  and  particle  'm'  at  the  contact  point  '  c' ,  as  shown 


schematically  in  Fig.  2.4,  are  given  by 

c  nm  in  .  n  .  , 4  m  m  .  n  n. 

A8  .  =  Au.  -  Au.  +  e.  Ao)  .  r.  -  A©  .  r, 
xii  ljk  j  k  j  k 

AOnm  m  n 
A0  .  =  Aw .  -  Aw  . 

ill 


(2.1) 

(2.2) 


where  the  quantity  e. =  the  permutation  symbols  used  in  tensor 

1]K 

representation  for  cross  product  of  vectors. 

The  relative  angular  rotation,  A0™,  is  related  to  the  contact 
couple.  The  rolling  and  torsional  resistance  at  contact  are  negligible 
for  elastic  particles  of  convex  shape,  but  significant  for  angular 
particles.  For  simplicity,  we  limit  our  discussion  to  particles  with 
convex  shape  and  neglect  the  effect  of  contact  couples. 

The  relative  displacement  representing  spring-stretch  at  the 
contact  point  'c',  A8™,  is  related  to  the  contact  force,  Af™,  by  a 
general  expression  as 

Af1™^  K™  A6 nm-  Af™  (2.3) 

i  i]  ]  i 
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Particle  m 


Fig.  2.4  Kinematics  of  two  rigid  particles  of  convex  shape 
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where  the  residual  force  Af™  and  the  contact  stiffness  tensor  K™  are 
described  below. 


2.3.1  Residual  Contact  Force 

During  a  small  strain  increment,  the  material  is  treated  to  be 
linear  governed  by  the  instantaneous  stiffness  tensor  K™.  As  a  result 
of  this  linearization,  the  resultant  contact  force  due  to  the  strain 
increment,  may  exceed  the  contact  strength.  The  effects  are  accounted 

A 

by  the  residual  force  Af . .  In  the  case  of  sliding,  the  residual  shear 

A 

force,  Af  ,  is  the  amount  of  contact  shear  force  in  excess  of  the 
contact  shear  strength.  In  the  case  of  particle  separation,  the 

A 

residual  normal  force,  Afn,  is  the  amount  of  contact  normal  force  in 
excess  of  the  contact  tensile  strength. 

A 

The  force  vector  Af^  for  this  contact  is  expressed  as  follows: 


Af.  =  Af  r.  +  Af  n. 
1  r  i  n  i 


(2.4) 


where  in  is  the  unit  vector  in  the  direction  of  resultant  shear  on  the 


contact  plane.  The  vector  in  is  normal  to  the  contact  plane. 


2.3.2  Contact  Stiffness  Tensor 

It  is  assumed  that  the  tangential  and  compressional  stiffness  are 
decoupled.  The  relationship  between  the  resultant  shear  force  f  and 
the  resultant  shear  displacement  8 r  on  the  contact  plane  is  fr=  5^. 
When  the  contact  force  exceeds  the  frictional  strength,  sliding  occurs 
and  vanishes  as  illustrated  in  Fig.  2.5.  The  relationship  between  the 
contact  force  f  and  the  relative  displacement  8^  in  the  normal 


A.  Normal 


Fig.  2.5  Force  and  displacement  relationships  in  the  directions  of:  a) 
normal,  and  b)  shear  on  a  plane  of  contact  between  two 
particles. 
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direction  is  f  =  k  5  .  When  the  contact  force  tends  to  be  in  tension, 
n  n  n 

particle  separation  occurs  and  k^  vanishes.  The  stiffness  tensor  at  the 
contact  between  particles  'n'  and  'm',  K™1,  takes  the  form 


K™  =  k  n.n.  +  k  (  ss  +  t.t.)  (2.5) 

where  n,  s  and  t  are  the  basic  unit  vectors  of  the  local  coordinate 
system  at  each  contact  as  shown  in  Fig.  2.6.  The  vector  n  is  normal  to 
the  contact  plane.  The  other  two  orthogonal  vectors  s  and  t  are  on  the 
contact  plane. 

A  frequently  used  formulation  for  contact  stiffness  between  two 
elastic  spheres  is  the  Hertz  theory  (Johnson  1985) .  Considering  the 
contact  area  to  be  circular  with  a  parabolic  pressure  distribution,  the 
deformation  at  the  contact  is  obtained  from  elasticity  solutions  for 
pressure  loads  on  semi-infinite  media.  This  leads  to  the  expression  of 
normal  stiffness  as  follows; 


kn  =  C1  rEp  <  2 


n  .  a 


(2.6) 


r  E 


y 

1/3 

where  is  a  constant  equal  to  3  ,  r  is  the  radius  of  particles,  a  = 

1/3,  E  =  G  /  ( 1 — v  ),  G  and  v  are  the  particle  shear  modulus  and 

y  y  y  y  y 

Poisson's  ratio  respectively . 

For  a  contact  of  two  topographically  smooth  spheres,  the  tangential 
stiffness  under  oscillating  contact  force  was  studied  by  Mindlin  and 
Deresiewicz  (1953) .  The  tangential  force  at  the  contact  results  in 
deformation  caused  by  the  development  of  slip  over  a  part  of  the  contact 
surface.  Clearly,  when  the  tangential  force  exceeds  the  frictional 
strength  at  the  contact,  sliding  takes  place.  A  general  expression  for 


Fig.  2.6  Local  co-ordinate  system  at  an  inter-particle  contact 
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the  tangential  stiffness  of  two  particles  can  be  written  as 
of  the  contact  force  and  the  particle  properties  as  follows: 


K  =  C_K  (1 
r  2  n 


where 


f  tan$ 
n 


a 


function 


(2.7) 


C  =  _ 

2  2-v 

P 

^  is  the  friction  angle  between  two  particles  and 
shear  force  at  the  contact. 


(2.8) 

is  the  resultant 


2.4  BEHAVIOR  AT  MICRO-ELEMENT  LEVEL 

Since  the  assembly  consists  of  a  large  number  of  particles,  it  is 
expedient  to  treat  translation  and  rotation  of  discrete  particles  as 
continuum  fields.  On  this  basis,  the  discrete  system  can  be  transformed 
into  an  equivalent  continuum  and  its  behavior  can  be  described  by  the 
continuum  concepts  of  stress  and  strain. 

The  strain  and  stress,  for  a  micro-element,  are  defined  in 
connection  with  the  relative  movement  of  the  particles  and  the  resulting 
contact  forces,  respectively.  The  stress-strain  relationship  is  then 
derived  based  on  contact  behavior. 


2.4.1  Strain 

At  the  micro-element  level,  movement  of  particles  can  be  determined 
by  displacement  gradient,  Au1^  . ,  with  the  usual  affine  (homogeneous 

1  r  1 

deformation)  assumption.  Similarly,  rotation  of  each  particle  can  be 
determined  by  the  gradient  of  particle  rotation.  However,  with  the 
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absence  of  couple  stress  of  the  micro-element,  the  gradient  of  particle 

rotation  can  be  neglected.  Therefore,  particles  within  this  group  are 

assumed  to  have  the  same  rotation,  A(on. 

1 

It  has  been  shown  (Chang  1989,  Chang  and  Ma  1990)  that,  for 
granular  material,  it  is  more  precise  to  define  an  asymmetric 
deformation  strain  in  terms  of  the  continuum  variables  u11  .  and  Ao>n  in 


the  following  manner: 

.  n  ,  n  i  -n 

Ac  .  .  *  Au.  .  +  e.  . .  A co. 
i]  i,]  i]k  k 


(2.9) 


The  symmetrical  part  of  Ac".  is  equal  to  the  symmetrical  part  of 


displacement  gradient,  i.e., 

.n  ,  n  1  . .  n  ,  ,  n  . 

U])  d/])  2  l,]  ],r 


(2.10) 


representing  the  usual  symmetric  Cauchy  strain  of  the  micro-element. 


The  skew  symmetric  part  of  Ac ^ .  is  given  by 


.  n  .  n  -n 

[i]]  [1,3]  31k  k 


(2.11) 


where  Au..  (Au.  Au11  . )  /  2  is  the  rigid  body  rotation  tensor.  The 

l  A,  3  J  Af]  ]/A 

skew  symmetric  part  of  Ac  ^  thus  represents  the  net  spin  of  particles 
(i.e.,  the  difference  between  rigid  body  rotation  of  the  micro-element 
and  the  average  rotation  of  particles) . 

The  strain  defined  in  this  way  not  only  is  theoretically  more 
generic  but  also  furnishes  a  convenient  kinematic  relationship  which 
relates  strain  to  the  spring  stretch,  AS™,  at  the  contact  between  the 


center  particle  '  n'  and  its  neighbor  particle  'm' ,  given  by 


,  t  nm  .  n  .  nm 
A6  ,  =  Ac  u  . 

A  x3  1 


(2.12) 


where,  L.  is  the  branch  vector  joining  the  centroids  of  particle  'n' 


.  nm  nc 


and  particle  'u' .  Note  that  L.  =  r." 

1  1 


me  ,  nc 

r .  ,  and  r .  =  the  vector 
3  ] 
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measured  from  the  centroid  of  the  particle  '  n'  to  the  contact  point  'o', 
as  illustrated  in  Fig.  2.4. 

It  is  noted  that  Eq.  2.12  is  based  on  a  kinematic  assumption  of 
affinity  (i.e.,  linear  displacement  field).  However,  the  displacement 
field  is  in  general  non-linear.  By  forcing  particles  to  move  in 
accordance  with  a  linear  displacement  field,  we  impose  a  constraint  to 
the  system.  As  a  result,  this  assumption  leads  to  a  stiffer  system. 

The  effect  is  particularly  pronounced  at  large  strain  level  when 
significant  amount  of  sliding  have  taken  place  (Chang  and  Misra  1990b) . 

In  this  chapter,  we  derive  a  non-affine  kinematic  relationship  so 
that  it  releases  the  restriction  due  to  linear  displacement  field 
assumption.  This  is  done  by  allowing  the  center  particle  to  move  freely 
to  achieve  equilibrium  of  contact  forces  while  the  surrounding  particles 
ars  fixed.  Thus,  in  addition  to  the  movement  in  accordance  with  the 
linear  displacement  field,  an  extra  movement,  is  expected  for  the  center 


particle.  This  movement  results  in  an  additional  stretch,  A6^,  at  all 
the  contacts  of  the  center  particle.  Thus  the  total  spring  stretch 


between  particles  '  n'  and  'm'  becomes: 

..nm  .  n  Tnm  ,  * 

M  .  =  Ae  .  .  L  .  +  Ad  . 

i  i]  ]  l 


(2.13) 


Since  the  center  particle  moves  to  satisfy  the  equilibrium  of 

* 

contact  forces,  the  extra  stretch  a6 .  can  be  obtained  from  the  force 


equilibrium  condition,  given  by 

„n  Ft  „nm:-i 
kn 


A6.  -  Tn.  Ae n 
l  ljm  jm 


T  =  -  F  K 

nm  \±  ic 


)  K™  L™ 

4*  v  m- 


(2.14) 


Eq.  2.14  shows  that  the  magnitude  of  A6  .  is  a  function  of  particle 
arrangement  and  contact  stiffness. 
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Substituting  this  equation  into  Eq.  2.13,  the  stretch  of  spring  at 


contact  relates  to  the  strain  of  the  micro-element  by 


AS™  =  (8  .  .L1™ 
l  lj  m 


Pn.  )  Ae n 
ljm  jm 


(2.15) 


When  the  displacement  field  is  linear,  T^_-m  vanishes  and  Eq.  2.15 
reduces  to  Eq.  2.12. 


2.4.2  Stress 

The  relationship  between  the  contact  forces  and  the  stress  of  the 
micro-element  can  be  defined  by  employing  the  theorem  of  stress  mean 
(Chang  and  Liao  1990) .  The  local  stress  Ao*^  is  expressed  as  the  volume 
average  of  the  dyadic  product  of  contact  force  Af™  and  branch  vector  L™ 

"iJ  1?  1  Li  lf:  l2'161 

where  Vn  is  the  volume  associated  with  the  n-th  micro-element. 

Summation  of  the  volume  Vn  over  all  micro-elements  is  equal  to  the  total 
volume  of  the  representative  unit,  such  that  V  =  ^  V. 


2.4.3  Local  constitutive  law 

Constitutive  equation  can  thus  be  established  to  relate  the  local 
strain,  Ae^,  and  the  local  stress,  Ac?j,  based  on  the  following 
relationships:  (1)  stress  versus  contact  forces  (Eq.  2.16),  (2)  force 
versus  spring-stretch  at  contact  (Eq.  2.3),  and  (3)  strain  versus  sprir.g- 
stretch  at  contact  (Eq.  2.15).  The  constitutive  equation  for  the  micro¬ 
element  can  thus  be  obtained  in  the  following  form: 

.  n  -  rn  .  n  AAn 
i]  ijkl  kl  lj 


(2.17) 
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where  C.  .  is  the  local  stiffness  tensor  for  the  n-th  micro-element 
i  lkl 


given  by 


cn. 


1  r  nm  „nm  .  c  . nm 

"  )  L»  •  xa  .  ( o  ,  Li 

n  L  i  “in  rV  i 


rn  ) 
"pkl 


(2.18) 


1]kl  2Vn  m  “]P  '”Fn 
The  instantaneous  local  stiffness  tensor  derived  in  Eq.  2.18  is  a  function 

of  contact  stiffness  and  the  position  of  particles.  The  constitutive 
tensor  has  the  following  properties  of  symmetry: 


,n 


,n 


p**  _  p“  _  pti  _  pn 

ijkl  jik.1  jilk~  ikjl‘ 

An 

The  stress  Ac^  in  Eq.  2.17  for  the  n-th  micro-element,  resulting 
from  the  residual  forces  due  to  sliding  and  separation,  is  given  as 
follows; 

Aon.  =  ^  (2.19) 

^  2Vnm  1  1 


2.5  BEHAVIOR  AT  REPRESENTATIVE  UNIT  LEVEL 

In  what  follows,  we  consider  a  given  volume  to  be  representative  of 

the  granular  solid  such  that  the  boundaries  of  the  said  volume  are 

subject  to  displacements  compatible  with  a  uniform  overall  strain. 

Under  such  conditions,  Hill  (1967)  has  shown  that,  for  heterogeneous 

material,  the  overall  stress  and  strain  can  be  expressed  as  the  volume 

averages  of  their  corresponding  quantities  at  local  level. 

Thus  the  overall  stress  and  strain,  denoted  by  Ao .  .  and  Ae  .  are 

ij  i]' 

regarded  as  volume  averages  of  the  local  stress  and  local  strain  at  the 
micro-element  level,  such  that 


Ao  .  . 

1 

vn 

n 

AO  .  . 

(2. 

.20) 

d] 

V 

n 

dj 

Ae 

1 

X 

vn 

n 

Ae 
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Corresponding  to  the  overall  soress  and  strain,  it  is  reasonable  to 
define  an  overall  stiffness  tensor  for  the  representative  unit  suoh  that 
the  overall  stress-strain  relationship  can  be  expressed  in  the  following 
form: 
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To  facilitate  a  rational  averaging  process,  we  introduce  a 

'distributive  law  of  strain'  which  relates  the  overall  strain  Ae  .  .  to 

i] 

the  strain  fraction  of  a  micro-element  Ae^,  given  by, 


.  n  _n  . 

Ae .  .  =  G.  .  Ae  ,  . 
lj  ijkl  kl 


(2.2" 


where  the  dimensionless  tensor  G.  .,  associated  with  each  micro- 

i  jkl 

element,  describes  indirectly  the  degree  of  heterogeneity  of  the 
material.  From  Eg.  2.27,  the  following  condition  must  be  satisfied: 


1  Gn. 
i  il 


jkl  Iijkl; 


(2.28) 


where  is  a  fourth  rank  identity  tensor  defined  in  terms  of 

Kronecker  delta  6 . .  as 
i] 

I.  ...-  i  (8  ..  8  ..+  8  ..8  ..  )  (2.29) 

i]kl  2  lk  ]1  ]1  ik' 

Thus,  mathematically,  i-s  a  probability  density/distribution 

function  in  tensorial  form. 

For  homogeneous  material,  the  strain  is  equally  distributed  to  each 


Therefore  the  distribution  tensor  G.  .  has  the 

i  jkl 


micro-element . 

following  form: 

.Ji 

Gn  . ,  .  =  I.  ...  ^  (2.30) 

ijkl  ijkl  V 

With  the  distributive  law  of  strain  (Eq.  2.27),  from  volume  average 
of  Eq.  2.17,  and  using  Eqs.  2.20  and  2.25,  the  stiffness  tensor  can  be 

written  as  a  'weighted'  volume  average  of  micro-element  stiffness  tensor 


Ci jki '  9iven  by 


.n 


C.  ...  -  )  C  ■  G  ,  . 

ijkl  i  ljmn  mnkl 


(2.31) 


Since  the  representac*ve  unit  consists  of  large  number  of  randomly 
arranged  particles,  the  heterogeneous  system  can  be  viewed  as  a 
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statistical  homogeneous  system.  The  strain  distribution  has  no 
correlation  with  the  geometric  location  and  the  distribution  tensor 
^ijkl  sk°u^  not  defined  in  geometric  space.  Instead,  strain 
distribution  is  largely  depending  on  the  distribution  of  stiffness 
^ijkl'  inc*icating  that  the  distribution  tensor  should  be  a  function 
defined  in  stiffness  space. 

It  is  thus  reasonable  to  establish  the  distribution  tensor  in  terms 

of  the  variance  of  the  micro-element  stiffness  Cn.,  ,  with  resort  to  the 

1  jkl 

average  stiffness  of  the  packing  defined  as  a  dimensionless 


tensor  given  by 

Vn  =  c~^  (Cn  -  C  ) 

i]kl  ljmn  mnkl  mnkl 


(2.32) 


The  distribution  tensor  can  be  related  to  in  many  forms.  It 

is  found  that  a  simple  and  suitable  form  for  the  distribution  tensor  is 


a  generalized  tensor  form  of  Maxwellian  type,  i.e., 

vn 

Gn  ..  ,=  A.  .  Exp  (-vn  .  . )  ^  (2.23) 

ljkl  ljmn  r  mnkl  V 

The  constant  tensor  A.  ..  .  involved  in  the  distribution  function  must 

i]kl 

satisfy  the  condition  of  identity  (i.e.,  Eq.  2.28).  Thus  it  can  be 
obtained  by 


iiki  -5  ^ExP<-vinki»vn 


where  the  exponential  term  can  be  expanded  into  a  series  as  follows: 


^  vmnkl^  ^mnkl  vmnkl+  2!  vra 


(2.35) 


r  mnkl  mnkl  mnkl  2!  mnpq  pqkl  . 

The  distributive  tensor  can  also  be  obtained  fron  a  self-consistent 

method  in  which  the  micro-element  is  considered  to  be  an  inhomogeneity 

with  constitutive  tensor  Cn.,.  embedded  in  an  infinite  medium  of 

ijkl 

constitutive  tensor  similar  to  that  used  in  the  study  of 


pclycrystalline  behavior  (Hutchingson,  1970,  Nemat -Nasser  and  Mahrabadi 
1984)  .  The  usual  self-consistent  method  can  not  be  directly  applied  to 
granular  materials  because  particles  in  granular  media  are  not  directly 
in  analogy  to  the  crystals  in  polycrystalline.  In  contrast  to 
polycrystalline  where  deformation  occurs  primarily  in  crystals,  the 
deformation  of  granular  material  principally  occurs  at  contacts  between 
particles.  In  addition,  the  fact  that  granular  material  has  an 
asymmetric  strain  tensor  adds  complexity  to  the  mathematical  process  in 
the  self-consistent  method.  (Misra,  1990) . 

2.6  STOMARY 

Perceiving  granular  material  as  a  collection  of  particles,  a 
constitutive  law  for  granular  material  is  derived  based  on 
micromechanics  approach,  taking  into  account  the  mechanisms  of  sliding 
and  separation  of  particles.  For  simplicity,  the  discussion  of  the 
present  theory  is  restricted  to  the  conditions  of  particles  with  convex 
shape . 

Summary  of  the  necessary  equations  for  defining  the  constitutive 
tensor  is  given  in  Table.  2.1.  Given  also  in  Table  2.1  are  the  require 
parameters  such  as  inter-particle  contact  properties  and  initial 
structure  of  micro-elements.  In  the  present  theory,  structure  of  the 
granular  material  is  intended  to  be  represented  statistically  by  a  set 
of  micro-elements  where  each  micro-element  is  defined  by  a  specific 
arrangement  of  particle  group.  This  chapter  is  focused  on  the 
micromechanics  frame  work  for  the  constitutive  law  of  granular  material 


CONSTITUTIVE  LAW: 


Representative  Unit  :  No.  of  Micro-Elements,  Volume  of 

Representative  Unit 


Each  Micro-Element  :  No.  of  Contacts,  Volume  of  Micro-Element 
Each  Contact  :  Branch  Vector,  Contact  Stiffness 


Table  2.1  :  Summary  of  the  constitutive  equations. 
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with  the  assumption  that  we  already  have  the  information  regarding  the 
representative  set  of  micro-elements. 

It  is  noted  that  the  constitutive  law  is  in  an  incremental  form. 

At  end  of  each  stress  increment,  the  initial  structure  of  these  micro¬ 
elements  are  updated.  Thus,  the  non-linearity  due  to  evolution  of 
microstructure  is  incorporated. 

The  unique  features  of  the  present  theory  are  listed  as  follows: 

1)  It  explicitly  accounts  for  the  effects  of  micro-structure,  thus  is 
capable  of  modelling  inherent  anisotropy  of  material  structure. 

2)  It  accounts  for  mechanical  properties  at  granular  contact  with  the 
consideration  of  mechanism  of  sliding  and  separation  between 
particles,  thus  is  capable  of  modelling  strength,  residual 
strength,  behavior  of  yielding,  and  behavior  of  plastic  flow  of 
the  material. 

3)  It  accounts  for  the  evolution  of  micro-structure,  thus  is 
capable  of  modelling  strain  dependency,  path  dependency,  stress- 
induced  anisotropy,  and  behavior  related  to  the  memory  of  the 
material. 

Therefore  the  present  theory,  in  a  more  fundamental  way,  models 


many  aspects  of  material  behavior.  This  constitutive  theory  is 
potentially  useful  for  comprehensive  modelling  of  complex  behavior  for 
granular  material. 
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CHAPTER  3 

PREDICTIONS  COMPARED  WITH  COMPUTER  SIMULATION 


The  present  micromechanics  model  is  illustrated  to  predict  the 
stress-strain  behavior  of  a  random  packing  of  planar  disks.  To  evaluate 
the  applicability  of  the  present  model,  the  predicted  results  are 
compared  with  that  obtained  from  computer  simulation  for  the  same  random 
packing. 

A  periodic  packing  is  used  in  this  example.  The  representative  unit 
of  the  packing  shown  in  Fig.  3.1  is  made  up  of  25  particles,  12 
particles  of  radius  0.21  mm  and  13  particles  of  radius  0.105  mm.  As 
already  noted,  each  particle  in  the  representative  unit  is  in  contact 
with  several  neighboring  particles.  These  particles  form  a  particle 
group  known  as  micro-element.  A  schematic  representation  of  the  25 
micro-elements  corresponding  to  each  particle  is  shown  in  Fig.  3.2.  Each 
micro-element  consists  of  a  reference  (or  central)  particle  in  contact 
with  its  neighbors.  Microstructure  of  the  random  packing  shown  in  Fig. 
3.1  is  represented  by  the  25  micro-elements. 

In  the  stress-strain  prediction,  the  contact  normal  stiffness  is 

taken  to  be  k  =200  kN/m.  The  range  of  the  stiffness  ratio,  k  /k  ,  is 
n  3  s  n 

from  0.3  to  1,  and  the  range  of  inter-particle  friction  angle  is  14 
to  25  degrees.  This  range  of  inter-particle  friction  angles  is 
comparable  to  the  measured  values  reported  for  dry  spheres  manufactured 
from  glass  or  steel  materials  (Skinner,  1969)  . 
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3.1  ONE  DIMENSIONAL  TEST 

The  stress-strain  behavior  of  the  packing  is  predicted  under  two 

loading  conditions:  one  dimensional  and  biaxial  compression.  In  the 

one-dimensional  compression  test,  the  packing  was  first  isotropically 

2 

compressed  to  o_=  203  kPa  (2.07  kg /cm  ).  .nen,  the  horizontal  strain  is 
held  constant,  i.e.,  At  =0,  while  the  samole  is  strained  in  the 

XX 

vertical  direction.  The  predicted  results  from  the  present  self- 
consistent  micromechanics  model  are  compared  with  that  obtained  from  the 
computer  simulation  method  (TRUBALL)  developed  by  Cundall  and  Strack 
(1979  a, b, c) . 

Fig.  3.3  shows  the  value  of  stress  ratio  Kq,  i.e.  c^/  o^.,  versus 

vertical  strain  for  various  values  of  From  the  state  of  isotropic 

compression,  Kq  decreases  from  1  sharply  to  a  value  and  stays  nearly 

constant  with  increasing  vertical  strain.  For  both  self-consistent  and 

computer  simulation  methods,  the  minimum  Kq  value  generally  increases 

with  decreasing  $  .  In  computer  simulation,  Kq  varies  from  0.27  to  0.16 

and  for  the  micromechanics  model,  it  varies  from  0.22  to  0.14  for  6  = 

U 

14°and  25°respectively .  Results  of  the  two  methods  show  reasonable 
agreement . 

3.2  BIAXIAL  COMPRESSION  TEST 

In  the  biaxial  loading  condition,  the  packing  was  first 

2 

isotropically  compressed  to  cq=  ^03  kRa  (2.07  kg/cm  ).  Then,  the 
horizontal  stress  is  held  constant  while  a  compressive  strain  is  applied 
in  the  vertical  direction.  The  results  of  the  prediction  are  plotted  in 
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term  of  stress  ratio,  q/p,  defined  as  follows:  q=(c  -  o, )  12,  and  d=  (o  + 

v  n  v 

o^)/2.  The  curves  of  stress  ratio  q/p  and  volumetric  strain  versus 

vertical  strain  are  shown  in  Fig.  3.4  for  $  =  25°  and  k  /k  =1.  Fiqure 

p  s  n 

3.5  plots  the  stress-strain  curves  for  k  /k  =0.3  and  $  =  25°. 

s  n  p 

Compared  with  Fig.  3.4,  the  effect  of  kg  on  strength  is  insignificant. 

The  predicted  stress-strain  curves  from  the  present  micromechanics 
model  show  reasonable  agreement  with  that  obtained  from  computer 
simulation  method.  The  discrepancies  are  caused  by  the  assumption  in 
the  micromechanics  model  that  the  micro-element  is  embedded  in  an 
equivalent  homogeneous  media  rather  than  in  a  heterogeneous  media.  The 
same  degree  of  discrepancies  was  observed  between  the  results  of  the 
micromechanics  model  and  the  computer  simulation  for  material  with 
larger  number  of  particles.  A  typical  example  predicted  using  the 
micromechanics  model  for  a  packing  with  250  particles  is  plotted  in  Fig. 

3.6  for  stress-strain  curves  at  three  confining  pressures.  Although  the 
example  is  for  idealized  packings  with  circular  disks,  it  is  noted  that 
the  stress-strain  response  resembles  very  much  the  behavior  of  sand. 

The  predicted  results  from  the  present  model  have  shown  reasonable 
agreement  with  that  obtained  from  computer  simulation  method.  The 
discrepancies  are  due  to  the  assumptions  inherent  in  the  self-consistent 
method  which  results  in  constraints  imposed  on  the  system.  However,  the 
general  behavior  remarkable  resembles  the  observed  behavior  for  sand. 

The  unique  features  of  the  present  theory  is  that  it  explicitly 
accounts  for  the  effects  of  micro-structure,  and  mechanism  of  sliding 
and  separation  between  particles.  Thus  the  model  is  capable  of 


Fig.  3.4  Stress-strain  curve  of  biaxial  test  for  k  /k  =  1  predicted  by 
computer  simulation  and  present  model.  s  n 


3 


45 


predicting  failure  strength  and  modelling  plastic  and  path-dependent 
stress-strain  behavior.  The  present  constitutive  theory  is  potentially 
useful  for  comprehensive  modelling  of  complex  behavior  for  granular 
material . 
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CHAPTER  4 

EVALUATION  WITH  EXPERIMENTAL  RESULTS 


Chapter  2  described  the  constitutive  model  for  granular  materials 
based  on  micro-scale  interaction  between  particles  including  the  effects 
of  sliding  and  separation  between  particles.  In  this  chapter,  the 
theory  is  preliminarily  evaluated  for  two  conditions,  namely,  small 
strain  and  large  strain. 

For  small  strain  conditions,  we  evaluate  the  capability  of  the 
constitutive  model  by  comparing  the  predicted  mechanical  properties, 
such  as  initial  moduli,  secant  moduli  and  damping  ratio,  with  those 
measured  from  experiments  on  sands  under  low  amplitude  cyclic  loading. 
The  focus  is  on  the  model  capability  in  simulating: 

1)  the  macroscopic  stress-strain  behavior  based  on  microscopic 
contact  property, 

2)  the  stress-induced  and  inherent  anisotropy, 

3)  the  response  and  energy  dissipation  of  material  under  cyclic 
loading,  and 

4)  the  effect  of  cemented  inter-particle  properties. 

For  large  strain  conditions,  the  predicted  stress-strain  behavior 
is  compared  with  experimental  results  on  different  types  of  sands  under 
various  loading  conditions,  such  as: 

1)  one  dimensional  compression, 

2)  cubical  tests, 

3)  circular  stress  paths  on  octahedral  plane,  and 
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4)  direct  simple  shear. 

The  comparison  of  salient  features  of  constitutive  behavior  is 
presented  to  evaluate  the  performance  of  the  developed  constitutive 
model  for  granular  material. 

4.1  MODEL  EVALUATION  FOR  SMALL  STRAIN  CONDITIONS 

The  assumption  of  uniform  strain  is  substantiated  with  the  results 
of  computer  simulation  fur  packings  under  small  strain  conditions  (Chang 
and  Misra  1989) .  With  this  assumption,  the  present  constitutive  model 
reduces  to  the  same  form  as  that  of  the  previously  developed  model 
(Chang  1988,  1989)  which  has  been  applied  to  the  studies  of  moduli  and 
damping  under  low  amplitude  cyclic  load.  A  few  examples  are  provided 
here  to  illustrate  the  model  capabilities. 

For  the  problems  of  small  strain,  the  structure  of  granular 
material  is  statistically  represented  by  a  set  of  binary-units  in 
various  orientations.  The  binary-unit,  defined  as  two  particles  in 
contact,  is  the  basic  constituent  unit  of  the  system.  For  convenience, 
in  all  examples  the  granular  material  are  envisioned  to  be  collections 
of  elastic  spheres.  Although  the  proposed  theory  is  not  limited  to 
spherical  particles,  it  is  difficult  to  implement  arbitrary  shaped 
particles  because  of  the  invalidity  cf  the  usual  Hertzian  contact  theory 
and  the  manifold  complexity  involved  in  the  description  of  packing 
geometry.  For  packings  with  equal-sized  spherical  particles,  the 
structure  can  be  characterized  by  the  distribution  of  contact 
orientations  and  the  average  coordination 


number . 
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4.1.1  Macro-behavior  in  terms  of  micro-property 

For  isotropic  packing  of  equal-sized  elastic  spheres,  the  shear 
modulus,  G,  and  the  Poisson's  ratio,  v,  were  derived  in  terms  of  contact 
properties  (Chang  et  al.  1990),  given  by 

V 


U 


10  -6v 


(4.1) 


G  =  F (z)  P (E  )  cr  (4.2) 

where  is  the  Poisson's  ratio  of  the  granule,  and  oc  is  the  confining 

stress.  In  Eq.  4.2,  F(z)  is  a  function  representing  packing  structure, 
given  by 

T  —  n 

(4.3) 

where  z  is  a  dimensionless  variable,  defined  by 
n 


F(z)  z1"0 


=  _M  3 

2  ~  3V  r  4jt  (1+e) 


(4.4) 


where  r  is  the  value  of  radii  of  spheres,  M/V  is  the  number  of  contacts 
(double  count)  per  unit  volume  of  the  representative  unit.  The 
parameter  z  is  also  a  function  of  the  void  ratio,  e,  and  the  average  co¬ 
ordination  number,  n.  The  void  ratio  e  is  defined  as  the  ratio  of  the 
volume  of  voids  to  the  volume  of  solids,  and  the  co-ordination  number  n 
is  defined  as  the  average  number  of  contacts  per  particle. 

In  Eq.  4.2,  P(E^)  is  a  function  of  granule  properties,  given  by 
P(E  )  =  C1(2+C2)E^"a  (4.5) 

In  general,  the  value  of  constant  Z^  is  between  0.5  to  1  and  the  value 
of  a  is  between  0.3  to  0.5.  For  Hertz-Mindlin  theory,  a  is  1/3,  C^= 

1.44,  E  =  G  /  ( 1— v  ),  G  and  v  are  the  shear  modulus  and  Poisson's  rati; 

P  u  p  p  P 

of  the  granule,  respectively.  C»  is  given  by 
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P 

Chang  et  al.  (1990)  have  shown  that  Eg.  4.2  becomes  identical  with 
Hardin  and  Black's  (1966)  expression,  using  a  =  0.5  and  an  appropriate 
relationship  of  e-n  curve.  It  has  also  been  shown  that,  for  a  =  0.33 
with  selected  values  of  modulus  and  Poisson's  ratio  for  quartz,  the 
predicted  results,  as  in  Fig.  4.1,  are  comparable  with  experimentally 
measured  shear  modulus  (Chung  et  al.  1984)  on  Monterey  No.  0  sand  in 
resonant  column  device  under  various  confining  pressures.  In  Fig.  4.1, 
the  predicted  results  are  also  compared  with  two  empirical  equations 
proposed  by  Hardin  and  Black  (1966),  and  Chung  et  al.  (1984). 

4.1.2  Stress-induced  and  inherent  anisotropy 

Chang  et  al.  (1990)  have  shown  that,  for  an  geometric  isotropic 
packing  under  isotropic  stress,  the  contact  forces  and  the  corresponding 
force-dependent  contact  stiffness  are  same  for  all  contacts.  As  a 
result,  the  constitutive  matrix  computed  from  the  contact  stiffness 
resembles  the  general  Hooke's  law  of  the  elastcity  theory  for  isotropic 
material . 

However,  when  the  stress  becomes  anisotropic,  the  contact  forces 
and  the  corresponding  force-dependent  contact  stiffness  are  no  longer 
same  for  all  contacts.  This  results  to  a  anisotropic  constitutive 
matrix  even  though  the  packing  geometry  is  isotropic.  The  change  of 
moduli  with  applied  stress  is  generally  referred  as  stress-induced 
anisotropy.  The  moduli  predicted  by  the  model  has  shown  good  agreement 


SHEAR  MODULUS  Gp  (pai  x  104) 
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CONFINING  PRESSURE  P  (pai) 


Fig.  4.1  Comparison  of  measured  and  predicted  results  of  shear  moduli  for 


sand  (Chang  et  al.  1991a) . 
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with  those  measured  from  experiments  conducted  by  Yanagisawa  (Chang  et 
al.  1990)  . 

Another  source  of  anisotropy  is  due  to  the  packing  anisotropy, 
termed  as  inherent  anisotropy.  Characterization  of  packing  structure  by 
fabric  tensor  was  discussed  by  Chang  and  Misra  (1990)  to  relate  material 
symmetries.  To  illustrate  the  model  capability  of  simulating  inherent 
and  induced  anisotropy,  predicted  moduli  are  compared  with  that  obtained 
experimentally  by  Stokoe  et  al.  (1985)  from  wave  velocities  measured  in 
various  directions  of  a  large  cubic  sample  of  sand. 

The  measured  initial  moduli  are  plotted  in  Fig.  4.2a  by  symbols, 
for  various  initial  confining  stress.  The  measured  moduli  are  different 
in  vertical  and  horizontal  directions,  showing  that  the  soil,  prepared 
by  air  pluviation  method,  has  an  inherent  transverse  isotropic  packing 
structure . 

After  the  sample  is  isotropically  stressed,  two  loading  conditions 
were  followed  to  show  the  stress-induced  anisotropy.  Figure  4.2b  shows 
the  results  of  increasing  vertical  stress  with  constant  horizontal 
stress.  Figure  4.2c  shows  the  results  of  decreasing  horizontal  stress 
with  constant  vertical  stress. 

The  predictions  in  all  cases  were  performed  (Chang  et  al.  1991a)  by 
using  the  material  properties,  G^,  and  and  a  parameter 
describing  the  distribution  of  contact  orientations.  The  comparison 
shows  that  the  model  can  effectively  simulate  both  inherent  and  induced 


anisotropy. 
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Fig.  4.2  Comparison  of  measured  and  predicted  shear  moduli  for  sand  under 

three  loading  conditions:  a)  isotropic  compression,  b)  oz 

increases,  0*0=  constant,  and  c)  0  *  constant,  0*0 
x  y  2  x  y 

decreases  (Chang  et  al.  1991a). 
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4.1.3  Cyclic  loading  and  related  energy  dissipation 

The  constitutive  model  is  capable  of  predicting  behavior  such  as, 
initial  and  secant  moduli,  and  damping  ratios,  under  low  amplitude 
cyclic  loading.  The  predicted  (Chang  et  al.  1990)  secant  moduli  and 
damping  ratio  are  compared  with  that  measured  by  Acar  and  El-Tahir 
(1986)  for  Monterey  No.  0  sand  and  shown  in  Figs.  4.3  and  4.4. 

In  the  Hertzian  contact  theory,  partial  sliding  can  occur  at  a 
contact  even  though  the  two  particles  are  still  intact  under  low 
amplitude  cyclic  loading.  The  partial  sliding  causes  degradation  of  the 
material  and  affects  the  value  of  secant  shear  modulus  Gs .  The  damping 
results  from  energy  dissipation  due  to  the  partial  sliding  between 
particles.  The  present  constitutive  theory,  with  microscopic  contact 
consideration,  is  particularly  useful  and  realistic  in  predicting  this 
type  of  behavior. 

4.1.4  Inter-particle  properties  -  cemented  sand 

The  model  is  capable  of  simulating  the  behavior  of  cemented  sand  by 
incorporating,  other  than  friction,  the  mechanism  of  adhesion  at  inter¬ 
particle  contact.  A  comparison  between  predicted  and  measured  results 
for  various  degrees  of  cementation  can  be  found  in  the  reference  of 
Chang  et  al.  (1990).  Typical  examples  are  shown  in  Figs.  4.3  and  4.4 
for  the  effects  of  cementation.  The  model  realistically  captures  the 
effects  of  cementation. 

4.2  MODEL  EVALUATION  FOR  FINITE  STRAIN  CONDITIONS 


MODULUS  RATIO  G„  /G 
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Fig.  4.3  Comparison  of  predicted  and  measured  secant  moduli  for  sa 


(Chang  et  al.  1990) . 
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4.4  Comparison  of  predicted  and  measured  damping  ratio  for  sand 
(Chang  et  al.  1990) . 
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For  the  problems  cf  small  strain,  the  structure  of  granular 
material  is  represented  statistically  by  a  set  of  binary-units  which  is 
not  sufficient  for  large  strain  problems.  Therefore,  the  structure  of 
granular  material  is  represented  statistically  by  a  set  of  micro¬ 
elements,  where  each  micro-element  is  defined  by  a  specific  arrangement 
of  particle  group. 

To  illustrate  the  current  model,  the  geometry  of  five  micro¬ 
elements  were  randomly  generated;  each  consists  of  eight  spherical 
particles  with  0.2mm  radius.  Each  micro-element  was  rotated  on  the  axes 
6  times  to  create  six  micro-elements  as  shown  in  Fig.  4.5,  resulting 
total  of  30  micro-elements.  The  idealized  material,  represented  by 
these  thirty  micro-elements,  has  material  symmetry  along  the  directions 
of  three  axes.  The  idealized  material  has  an  average  co-ordination 
number  of  7  and  void  ratio  0.7,  representing  an  uniformly  graded  medium 
dense  sand  with  rounded  particles. 

To  evaluate  the  capability  of  the  constitutive  model  at  large 
strain  conditions,  the  stress-strain  behavior  of  this  idealized  material 
is  predicted  for  the  following  loading  conditions: 

1)  one  dimensional  compression, 

2)  cubical  tests, 

3)  circular  stress  paths  on  octahedral  plane,  and 

4)  direct  simple  shear. 

The  predicted  results  are  compared  with  the  measured  behavior  from 
experiments  with  various  types  of  sands.  Since  the  predictions  for  all 
cases  are  based  on  the  idealized  material,  the  predictions  are  not 
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intended  to  fit  any  particular  experiment  curve.  The  compariscn  is 

merely  for  the  purpose  of  evaluating,  in  a  qualitative  sense,  the 

general  model  performance  under  various  loading  conditions. 

Only  three  basic  material  parameters  are  used  for  this  model, 

namely,  the  normal  contact  stiffness  k^,  the  shear  contact  stiffness  k  , 

and  the  inter-particle  friction  angle  $  .  For  the  predictions  of  all 

cases  listed  above,  the  normal  and  shear  contact  stiffness,  k  =  k  =  105 

n  s 

N/m,  are  used.  Although  the  non-linear  Hertzian  type  of  contact  can  be 
readily  used  for  these  predictions,  constant  contact  stiffness  are 
purposely  used  here  to  demonstrate  that  the  nonlinear  stress-strain 
behavior  is  predominantly  caused  by  the  mechanism  of  particle  sliding. 

4.2.1  One  dimensional  compression 

The  typical  behavior  of  granular  material  in  one  dimensional 
compression  is  shown  in  Fig.  4.6.  The  value  of  Kq  is  nearly  constant 
during  loading  but  gradually  increases  during  unloading.  At  microscopic 
level,  a  substantial  amount  of  sliding  occurs  during  loading  which 
results  locked-in  stress  during  unloading,  thus  reflects  by  a  higher  K  . 
The  present  theory,  accounting  for  the  mechanism  of  sliding  at  contact, 
is  capable  of  predicting  the  familiar  phenomenon  of  Kq  increase  with  the 
over-consolidation  ratio  (OCR) . 

The  idealized  material  is  also  used  to  study  the  one-dimensional 

compression  behavior  of  granular  material  preceded  by  isotropic 

0 

compression.  For  this  purpose,  $.,=  IT  is  selected  and  the  prediction 
is  compared  with  the  experimental  results  on  medium  dense  Napa  Basalt 
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(Lade  19"5)  shown  in  Fig.  4.7.  The  prediction  gives  excellent  agreement 
with  the  experimental  behavior  even  though  the  soil  structure  is 
hypothetical . 

The  behavior  of  granular  material  under  one-dimensional  compression 
preceded  by  anisotropic  consolidation  is  also  studied.  For  convenience, 
$  =  15  is  selected  for  the  predictions  of  three  initial  anisotropic 
stress  ratio,  0.5,  1.0  and  1.5.  The  results  are  in  good  agreement  with 
the  experimental  results  on  Tovoura  sand  (Okochi  and  Tatsuoka  1984),  as 
shown  in  Fig.  4.8.  It  may  be  seen  from  Fig.  4.8  that  the  Kq  values 
measured  under  different  initial  stress  conditions  converge  to  an 
identical  value  at  large  axial  stress. 


4.2.2  Cubical  tests 

The  constitutive  model  is  used  to  predict  stress-strain  and  volume 
change  behavior  of  the  idealized  granular  material  under  cubical 
triaxial  loading  conditions  with  different  values  of  b,  where  b  is 
defined  as 


<t„—  c. 


b  = 


V  03 


(4.7) 


Figure  4.9  shows  the  predicted  stress-strain  and  volume  change  behavior 
using  $  equal  to  25°.  The  predicted  behavior  is  compared  with  the 
experimental  results  obtained  from  cubical  triaxial  tests  on  dense 
Monterey  No.  0  sand  (Lade  and  Duncan  1973),  shown  in  Fig.  4.10.  It  is 
noted  that  the  nonlinear  stress-strain  behavior  is  predominantly  caused 
by  sliding  cf  particles  since  the  contact  stiffness  are  constants. 


6 


Fig.  4.7  Comparison  of  predicted  Kq  behavior  with  the  experimental 

results  on  medium  dense  Napa  Basalt  performed  by  Lade  (1975) 


Vertical  Stress  (kPa) 

Fig.  4.8  Comparison  of  predicted  Kq  behavior  with  the  experimental 

results  on  Toyoura  sand  for  different  initial  anisotropic  stress 
conditions  performed  by  Okochi  and  Tatsuoka  (1984)  . 


a3)  kg/cm  (  a1-a3)  kg/cm  k9/ 
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Although  the  prediction  is  based  or.  the  idealized  material 
represented  by  30  micro-elements,  the  predicted  trend  of  stress-strain 
and  volumetric  strain  characteristics  are  remarkably  similar  to  that  of 
experimental  values  for  Monterey  sand. 

The  predicted  variation  of  peak  frictional  angle  with  b  is  shown  in 
Fig.  4.11,  along  with  other  experimental  measurements  obtained  from 
literature  for  various  types  of  granular  materials  (Reades  and  Green 
1974,  Sutherland  and  Mesdarv  1969,  Lade  and  Duncan  1973) .  The  present 
theory,  even  with  the  idealized  material,  predicts  reasonable  results 
which,  both  qualitatively  and  quantitatively,  fall  in  the  range  of  that 
observed  from  experiments. 

The  predicted  failure  surface  on  octahedral  plane  is  compared  with 
the  failure  surfaces  empirically  hypothesized  by  Mohr-Coulomb,  Matsuoka 
(1974)  and  Lade  and  Duncan  (1973),  as  shown  in  Fig.  4.12.  It  is  the 
unique  capability  of  the  present  theory  which  predicts  the  failure 
surface  of  granular  material  in  a  stress  space  solely  on  the  basis  of 
inter-particle  friction  property. 

The  predicted  directions  of  the  strain  increments  form  acute  angles 
with  the  failure  surface  on  tnaxial  plane,  and  nearly  perpendicular  to 
the  failure  surface  on  octahedral  plane,  as  shown  in  Fig.  4.13.  This 
behavior  is  also  in  agreement  with  that  observed  experimentally  for  sand 
(Ko  and  Scott  1967,  Lade  and  Duncan  1973) . 


4.2.3  Circular  stress  paths 
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4.11  Predicted  variation  of  peak  friction  angle  with  b  for  the 


idealized  material,  along  with  other  experimental  measurements 
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ig.  4.12  Predicted  failure  surface  on  octahedral  plane  compared  with  the 
failure  surfaces  empirically  hypothesized  by  Mohr-Coulomb, 
Matsuoka  (1974)  and  Lade  and  Duncan  (1973) . 


i  directions  of  the  strain  increment  vectors  at  failure 
Laxial  plane  and  b)  octahedral  plane  for  the  idealized 
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To  further  evaluate  the  capability  of  the  model,  the  idealized 
material  subjected  to  circular  stress  path  on  octahedral  plane  is 
predicted.  A  classical  model,  such  as  an  isotropic  hardening  plasticity 
model,  would  predict  elastic  response  for  the  circular  stress  path 
within  the  yield  surface.  This  is  in  contradiction  to  the  observed 
experimental  results.  The  present  theory  models  the  plastic  sliding  at 
contacts  instead  of  using  a  priori  empirically  hypothesized  yielding 
surface.  Thus  the  present  model  is  capable  of  predicting  material 
response  under  complex  loading  conditions  such  as  the  circular  stress 
path.  Figure  4.14  shows  the  comparison  of  predicted  strain  path  with 
that  measured  from  experiments  on  Hostun  sand  performed  by  Lanier  and 
Zitouni  (1988) . 

4.2.4  Direct  simple  shear 

Predictions  of  strain-controlled  direct  simple  shear  are  performed 
on  the  idealized  material  using  equal  to  25°.  The  packing  is 
initially  consolidated  one-dimensionally  to  a  vertical  stress  of  134.4 
kPa,  followed  by  a  shear  strain.  Figure  4.15a  shows  the  predicted 
stress-strain  curves.  The  predicted  stress  ratio  at  failure  is  0.78. 

Figure  4.15b  shows  that  the  predicted  intermediate  principal  stress 
is  not  equal  to  the  mean  of  the  other  two  principal  stresses,  opposed  t: 
Hill's  (1950)  proposition  for  ideal  plastic  material.  However  the 
predicted  trends  are  in  agreement  with  that  measured  from  simple  shear 
tests  on  sand  reported  by  Resect  et  al.  (1967) . 


/0.5(o,1  +0- 
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simple  shear  test;  b)  Variati 
other  two  principal  stresses. 
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To  evaluate  the  behavior  of  rotation  of  principal  stress,  the 
predicted  results  are  compared  with  an  empirical  equation  by  Oda  and 
Konishi  (1974) .  The  empirical  equation  relates  the  inclination  angle 
y,  i.e.  angle  between  the  horizontal  plane  and  the  plane  on  which  the 
major  principal  stress  acts,  to  the  stress  ratio  in  simple  shear  test  on 
the  granular  material,  given  by 

t /g  =  k  tan  y  (4.8) 

This  equation  has  been  validated  for  Leighton  Buzzard  sand  in  monotonic 
and  cyclic  simple  shear  tests.  The  measured  values  of  k  for  Leighton 
Buzzard  sand  varies  from  0.58  to  0.67  (Cole  1967,  Budhu  1979).  Figure 
4.16a  shows  the  predicted  relationship  between  the  stress  ratio  and  tan 
V,  which  gives  a  value  of  k  =  0.64. 

In  addition  to  the  angle  v,  the  following  angles  p,  <o,  *,  and  \  are 
defined  respectively  as  the  angles  between  the  horizontal  plane  and  the 
plane  of  maximum  shear  stress,  the  plane  of  maximum  stress  ratio,  the 
plane  of  major  principal  stress  increment,  and  the  plane  of  the  major 
principal  strain  increment. 

The  predicted  results  of  these  angles  are  plotted  in  Fig.  4.16b 
which  shows  that  : 

a)  The  angle  \  is  not  equal  to  the  angle  v  indicating  that  the 
principal  axes  of  the  strain-increment  dc  not  coincide  with  the 
principal  axes  of  stress.  Thus  behavior  of  sand  is  in  contrary  to 
the  postulate  in  plasticity  theory  (Hill  1 9 ^ 0 ) . 

b)  The  angle  \  is  not  equal  to  the  angle  \  indicating  that  the  axes  of 
stress-increment  (%)  and  strain-increment  (i )  dc  not  coincide, 
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Fig.  4.16  a)  Predicted  relationship  between  the  stress  ratio  and  tan  v;  b) 

Predicted  results  of  angles  defined  as  the  inclinations  of  the  planes 
of:  1)  major  principal  stress,  y,  2)  maximum  shear  stress,  p,  3) 
maximum  stress  ratio,  co,  4)  major  principal  strain  increment,  and 
5)  major  principal  stress  increment,  x- 


which  is  the  characteristics  for  eiastic  material  (Hill,  1950) . 

Thus  sand  : s  not  an  elastic  material,  even  in  the  early  stage  of 
the  test. 

c)  The  angle  to  decreases  as  shear  strain  increases  and  approaches  to 
zero  at  large  shear  strain.  This  indicates  that  the  plane  of 
maximum  stress  ratio  does  not  occur  on  the  horizontal  plane  except 
at  a  later  stage  of  the  test. 

d)  The  angle  {5  is  approximately  zero  at  the  peak  shear  stress, 
indicating  that  the  horizontal  plane  is  the  plane  of  maximum  shear 
stress.  On  this  basis  the  peak  friction  angle  of  the  idealized 
material  in  direct  simple  shear  is  48.6°.  On  the  other  hand,  by- 
treating  the  horizontal  plane  as  the  maximum  stress  ratio,  the 
angle  of  friction  becomes  3 7C.  in  cubical  tests  the  maximum  value 
of  peak  friction  angle  is  48c  at  b  =  0.6,  whereas  for  triaxial 
compression  and  extension  cases  the  friction  angles  are  41°  and  45” 
respective.]../. 

The  predicted  trends  of  these  angles  are  in  agreement  with  that  measured 
from  simple  shear  tests  on  sand  reported  by  Rcscoe  et  al.  (1967) . 

4.3  SUMMARY 

The  developed  constitutive  model  has  been  evaluated  by  comparing 
the  predicted  behavior  of  the  idealized  material  with  that  obtained  from 
experiments  on  different  sands  under  various  loading  conditions.  For 
small  strain  conditions,  the  model  is  evaluated  for  its  applicability  ir. 
predicting  initial  moduli,  secant  moduli  and  damping  ratio  under  low 
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amplitude  leading, 
for  its  capability 


For  large  strain  conditions,  the  model  is  evaluated 
in  predicting  stress-strain-strength  behavior  under 


various  stress  paths. 

The  results  show  that  the  model  has  the  following  capabilities: 

1)  The  macroscopic  stress-strain  behavior,  such  as  moduli  and  failure 
strength,  can  be  predicted  in  terms  of  microscopic  contact 
properties . 

2)  The  effect  of  microstructure,  such  as  inherent  anisotropy,  can  be 
explicitly  accounted. 

3)  The  model  can  account  for  the  behavior  associated  with  particle 
interaction  and  sliding  at  contact,  such  as  stress-induced 
anisotropy,  path  dependency,  locked-in  stress,  plastic  flow, 
dilatancy,  friction  loss  (damping),  and  non-coaxial  behavior  under 
rotation  of  principal  stress. 

4)  The  model  can  account  for  materials  with  cemented  inter-particle 
properties . 

It  is  noted  that  in  all  the  predictions,  only  three  basic  material 
parameters  are  used  to  represent  the  stiffness  and  inter-particle 
friction  of  the  contact. 

The  theory  used  in  this  chapter  focuses  on  the  development  of 
mechanics  of  granular  material.  It  does  not  address  issues  on 
representation  of  geometric  structure  of  granular  material.  However, 
the  predicted  behavior  for  an  idealized  material  represented  by  30 
randomly  generated  micrc-elemer.cs,  has  shown  remarkable  similarity  to 
the  stress-strain  behavior  observed  from  experiments  under  various 
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loading  conditions.  This  agreement  in  the  predicted  and  measured 
behavior  indicates  the  potential  applicability  of  this  model.  Further 
study  is  under  progress  on  statistical  representation  of  microstructure 
of  granular  soil. 

The  purpose  of  this  chapter  is  to  illustrate,  aside  from  the 
prevailing  models,  an  alternative  way  of  modelling  stress-strain 
behavior  for  granular  material.  The  proposed  theory  is  conceptually 
simple  and  is  useful  in  discerning  the  essence  of  all  entities 
comprising  the  behavior  of  granular  material. 
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CHAPTER  5 

SUMMARY  AND  CONCLUSION 


Perceiving  granular  material  as  a  collection  of  particles,  a 
constitutive  law  for  granular  material  is  derived  based  on 
micromechanics  approach,  taking  into  account  the  mechanisms  of  sliding 
and  separation  of  particles.  In  the  present  theory,  structure  of  the 
granular  material  is  intended  to  be  represented  statistically  by  a  set 
of  micro-elements  where  each  micro-element  is  defined  by  a  specific 
arrangement  of  particle  group. 

The  constitutive  law  is  in  an  incremental  form.  At  the  end  of  each 
stress  increment,  the  initial  structure  of  these  micro-elements  are 
updated.  Thus,  the  non-linearity  due  to  evolution  of  micrcstructure  is 
incorporated. 

The  unique  features  of  the  present  theory  are  listed  as  follows: 

1)  It  explicitly  accounts  for  the  effects  of  micro-structure,  thus 
is  capable  of  modelling  inherent  anisotropy  of  material 
structure . 

2)  It  accounts  for  mechanical  properties  at  granular  contact  with 
the  consideration  of  mechanism  of  sliding  and  separation  between 
particles,  thus  it  capable  of  modelling  strength,  residual 
strength,  behavior  of  yielding,  and  behavior  of  plastic  flow  of 


the  material. 
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3)  It  accounts  for  the  evolution  of  micro-structure,  thus  is 

capable  of  modelling  strain  dependency,  path  dependency,  stress- 
induced  anisotropy,  and  behavior  related  to  the  memory  of  the 
material . 

Therefore  the  present  theory,  in  a  more  fundamental  way,  models  many 
aspects  of  material  behavior. 

The  developed  constitutive  model  has  been  evaluated  by  comparing 
the  predicted  behavior  of  the  idealized  material  with  that  obtained  from 
experiments  on  different  sands  under  various  loading  conditions.  For 
small  strain  conditions,  the  model  is  evaluated  for  its  applicability  in 
predicting  initial  moduli,  secant  moduli  and  damping  ratio  under  low 
amplitude  loading.  For  large  strain  conditions,  the  model  is  evaluated 
for  its  capability  in  predicting  stress-strain-strength  behavior  under 
various  stress  paths. 

The  results  show  that  the  model  has  the  following  capabilities: 

1)  The  macroscopic  stress-strain  behavior,  such  as  moduli  and  failure 
strength,  can  be  predicted  in  terms  of  microscopic  contact 
properties . 

2)  The  effect  of  microstructure,  such  as  inherent  anisotropy,  can  be 
explicitly  accounted. 

3)  The  model  can  account  for  the  behavior  associated  with  particle 
interaction  and  sliding  at  contact,  such  as  stress-induced 
anisotropy,  path  dependency,  locked-in  stress,  plastic  flow, 
dilatancy,  friction  loss  (damping) ,  and  non-coaxial  behavior  under 
rotation  of  principal  stress. 
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4)  The  model  can  account  for  materials  with  cemented  inter-particle 
properties . 

The  predicted  behavior  for  an  idealized  material  has  shown 
remarkable  similarity  to  the  stress-strain  behavior  observed  from 
experiments  under  various  loading  conditions.  This  agreement  in  the 
predicted  and  measured  behavior  indicates  the  potential  applicability  of 
this  model.  Further  study  is  under  progress  on  statistical 
representation  of  microstructure  of  granular  soil. 

The  results  illustrate  that  the  micromechanics  approach  of 
modelling  stress-strain  behavior  for  granular  material  is  conceptually 
simple  and  is  useful  in  discerning  the  essence  of  all  entities 
comprising  the  behavior  of  granular  material. 
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